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Abstract 

ON \ In this article, we assume that there exists a scalar D*D* molecular state in the 

J/i/xf) invariant mass distribution, and study its mass using the QCD sum rules. The 
predictions depend heavily on the two criteria (pole dominance and convergence of 
the operator product expansion) of the QCD sum rules. The value of the mass is 
about M D ,jj, = (4.43 ± 0.16) GeV, which is inconsistent with the experimental data. 
The D*D* is probably a virtual state and not related to the meson F(4140). Other 
ON . possibility, such as a hybrid charmonium is not excluded. 

^ PACS number: 12.39.Mk, 12.38.Lg 

Q_i' Key words: Molecular state, QCD sum rules 

1 Introduction 

^-j- \ Recently the CDF Collaboration observed a narrow structure (Y(4140)) near the J/ip(f> 

threshold with statistical significance in excess of 3.8 standard deviations in exclusive 
B + — > J/-00^ + decays produced in pp collisions at yfs = 1.96 TeV The mass and 
width of the structure are measured to be 4143.0 ± 2.9 ± 1.2 MeV and 11.7±|;o ± 3.7 MeV, 
respectively. The meson K(4140) is very similar to the charmonium-like state Y(3930) near 
the J/tpuj threshold [IE]. The mass and width of the Y(3930) are 3914.6+|;| ± 2.0 MeV 
^ ! and 34+f ± 5 MeV, respectively [3J. 

In Ref.jl], Liu et al study the narrow structure Y(4140) with the meson-exchange 
model, and draw the conclusion that the Y(4140) is probably a D*D* molecular state 
with J PC = ++ or 2++ while the V(3930) is its D*D* molecular partner. In Ref.[5], 
Mahajan argues that it is likely to be a D*D* molecular state or an exotic 
C$ ' hybrid charmonium. 

The mass is a fundamental parameter in describing a hadron, in order to identify 
the y(4140) as a scalar molecular state, we must prove that its mass lies in the region 
(4.1 — 4.2) GeV. In this article, we assume that there exists a scalar D*D* molecular state 
in the J/ip(ft invariant mass distribution, and study its mass with the QCD sum rules [6l [7]. 

In the QCD sum rules, the operator product expansion is used to expand the time- 
ordered currents into a series of quark and gluon condensates which parameterize the 
long distance properties of the QCD vacuum. Based on the quark-hadron duality, we can 
obtain copious information about the hadronic parameters at the phenomenological side 

The article is arranged as follows: we derive the QCD sum rules for the mass of the 
V(4140) in section 2; in section 3, numerical results and discussions; section 4 is reserved 
for conclusion. 
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2 QCD sum rules for the molecular state F(4140) 



In the following, we write down the two-point correlation function H(p) in the QCD sum 
rules, 

n(p) = if d*xj p - x (0\T{j(x)j\0)} |0) , (1) 
J(x) = c{x)-f^s{x)s(x)^c(x) , (2) 

we choose the scalar current J(x) to interpolate the molecular state Y(4140). 

We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operator J(x) into the correlation function H(p) to obtain the 
hadronic representation [HE]. After isolating the ground state contribution from the pole 
term of the Y(4140), we get the following result, 

n « ■■ (3) 

where the pole residue (or coupling) Ay is defined by 

x Y = (o|j(o)|y( P )). (4) 

In the following, we briefly outline the operator product expansion for the correlation 
function H(p) in perturbative QCD. The calculations are performed at the large space-like 
momentum region p 2 <C 0. We write down the "full" propagators Sij(x) and Cij(x) of a 
massive quark in the presence of the vacuum condensates firstly [7], 

+ l -^ m s(sg s aGs) * - ^^G%{^ + a*"*) + • - - , (5) 

r f v j f^^-ik xj &ij 9sG?f a a p()6 + m c ) + Qt + m c )<T af} 

) ~ (2vr)4 J d Ke \)t-m c 4 (k 2 - m 2 ) 2 

ir 2 a s GG k 2 + m c }6 \ 

+ T { — )6 ^ mc (k 2 -m 2 r + "7' (6) 

where (sg s aGs) = (sg s a a pG a ^ s) and ( as !? G ) = ( aaGa P G — ^ then contract the quark fields 
in the correlation function H(p) with Wick theorem, and obtain the result: 

n(p) = i f d^xe^Tr^^x^C^-x^TrirC^ix^Sn^-x)] , (7) 

where the i, j, m and n are color indexes. 

Substitute the full s and c quark propagators into the correlation function II (p) and 
complete the integral in the coordinate space, then integrate over the variables in the 
momentum space, we can obtain the correlation function H(p) at the level of the quark- 
gluon degrees of freedom. 
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We carry out the operator product expansion to the vacuum condensates adding up to 
dimension-10 and take the assumption of vacuum saturation for the high dimension vac- 
uum condensates, they are always factorized to lower condensates with vacuum saturation 
in the QCD sum rules, and factorization works well in large N c limit. In calculation, we 
observe that the contributions from the gluon condensate are suppressed by large denom- 
inators and would not play any significant roles [H El [TOl E] • 

Once analytical results are obtained, then we can take the quark-hadron duality and 
perform Borel transform with respect to the variable P 2 = —p 2 , finally we obtain the 
following sum rule: 



so 



4(m c +m s ) 2 



dsp(s)e 'm 7 , 



(8) 



where 



p(s) = Po(s) +p {Ss) (s) + pf GG) (s) + pf GG) (s) (^—) + p { ss)'4s), (9) 



the lengthy expressions of the spectral densities po(s), P(g S )(s), Pf GG )( s )' P(GG)( S ) an< ^ 
P{ss) 2 ( s ) are presented in the appendix. 

Differentiating the Eq.(8) with respect to -X^, then eliminate the pole residue Ay, we 
can obtain a sum rule for the mass of the y(4140), 



My 



4K +ms F & d(-l/Mf ( S ) e g 



(10) 



3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = —(0.24 ± 0.01 GeV) 3 , 
(ss) = (0.8±0.2)(gg), (sg s aGs) = rog(sa), = (0.8 ± 0.2) GeV 2 , = (0.33 GeV) 4 , 

m s = (0.14 ± 0.01) GeV and m c = (1.35 ± 0.10) GeV at the energy scale p = 1 GeV 

In the conventional QCD sum rules [6j[7], there are two criteria (pole dominance and 
convergence of the operator product expansion) for choosing the Borel parameter M 2 and 
threshold parameter so- 

In Fig.l, we plot the contribution from the pole term with variation of the threshold 
parameter sq. From the figure, we can see that the value sq < 20 GeV 2 is too small to 
satisfy the pole dominance condition. 

In Fig. 2, we plot the contributions from different terms in the operator product ex- 
pansion. The contribution from the term ) is very small, the contributions from 
the terms involving the gluon condensates are less than (or equal) 10% at the values 
M 2 > 2.3 GeV 2 and s > 22 GeV 2 , the gluon condensate plays a minor important role. 
The vacuum condensates of the highest dimension (ss}(sg s aGs)+(sg s (jGs) 2 serve as a 
criterion for choosing the Borel parameter M 2 and threshold parameter so- At the val- 
ues M^ in > 2.6 GeV 2 and sq > 23 GeV 2 , their contributions are less than 15%. The 
contribution from the vacuum condensate of high dimension (ss) 2 varies with the Borel 
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Figure 1: The contribution from the pole term with variation of the Borel parameter M 2 . 
The notations a, (3, 7, A, r, £ and p correspond to the threshold parameters sq = 19 GeV 2 , 
20GeV 2 , 21 GeV 2 , 22 GeV 2 , 23 GeV 2 , 24 GeV 2 and 25 GeV 2 , respectively. 

parameter M 2 remarkably and serves as another criterion for choosing the Borel parame- 
ter M 2 and threshold parameter s . At the value s > 23 GeV 2 and M 2 > 2.4 GeV 2 , its 
contribution is less than (or equal) 20%. In the region M 2 > 2.6 GeV 2 , the leading con- 
tribution comes from the perturbative term, while the next-to-leading contributions come 
from the terms (ss) + (sg s aGs). The operator product expansion is convergent at the val- 
ues M^ in > 2.6 GeV 2 and sq > 23 GeV 2 . For the central values of the input parameters, 
the contribution from the pole term is larger than 49% at the values M^ ax < 3.0 GeV 2 
and s > 23 GeV 2 . 

In this article, the threshold parameter and the Borel parameter are taken as so = 
(24 ± 1) GeV 2 and M 2 = (2.6 - 3.0) GeV 2 respectively, the contribution from the pole 
term is about (49 — 72)% for the central values of the other input parameters, the two 
criteria of the QCD sum rules are full filled [61 [7J. One may expect to take smaller Borel 
parameter and threshold parameter to satisfy the two criteria of the QCD sum rules 
marginally, however, it is not feasible. The contributions from the different terms in the 
operator product expansion change quickly with variation of the Borel parameter M 2 at 
the value M 2 < 2.6 GeV 2 (see Fig. 2) and will not result in a stable sum rule for the mass 
(see Fig.3). 

Taking into account all uncertainties of the input parameters, finally we obtain the 
values of the mass and pole residue of the Y, which are shown in Figs. 4-5, 

M Y = (4.43 ± 0.16) GeV , 

Ay = (5.46 ± 1.21) x 10~ 2 GeV 5 . (11) 

The central value My = 4.43 GeV is about 200 MeV above the D*D* threshold, the D*D* 
is probably a virtual state and not related to the meson Y(4140). Other possibility, such 
as a hybrid charmonium is not excluded. We can explore the hidden charm two-body 
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Figure 2: The contributions from the different terms with variation of the Borel pa- 
rameter M 2 in the operator product expansion. The A, B, C, D, E and F correspond 
to the contributions from the perturbative term, (ss) + {sg s aGs) term, ( as ) term, 
(^§^-)+{^§^-) [(ss) + (sg s aGs) + (ss) 2 ] term, (ss) 2 term and (ss)(sg s aGs)+ (sg s aGs) 2 
term, respectively. The notations a, (3, 7, A, r and p correspond to the threshold parame- 
ters s = 20 GeV 2 , 21 GeV 2 , 22 GeV 2 , 23 GeV 2 , 24 GeV 2 and 25 GeV 2 , respectively. Here 
we take the central values of the input parameters. 
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decay J/ip(j) an d the open charm two-body decays D S D S , D S D* to make further studies. 
More experimental data are still needed. 

From Eq.(ll), we can see that the uncertainty of the mass My is rather small (about 
3.6%) while the uncertainty of the pole residue Ay is rather large (about 22.2%). The 
uncertainties of the input parameters ({qq), {ss), {sg s aGs), m s and m c ) vary in the range 
(7 — 25)%, so the uncertainty of the pole residue Ay is reasonable. We obtain the value 
of the mass My through a fraction (see Eq.(lO)), the uncertainties in the numerator and 
denominator which origin from a given input parameter (for example, {ss), (sg s crGs)) 
cancel out with each other. It is not unexpected that the net uncertainty is smaller than 
the uncertainties of the input parameters. 

At the energy scale /i = lGeV, ^ « 0.19 p3], if the perturbative 0{a s ) corrections 
to the perturbative term are companied with large numerical factors, 1 + £(s,m c )^-, for 
example, £(s, m c ) > — 5, the contributions may be large. We can make a crude estima- 
tion by multiplying the perturbative term with a numerical factor, say 1 + £(s, m c )^- = 2, 
the mass My decreases slightly, about 20MeV, the pole residue Ay increases remark- 
ably. The main contribution comes from the perturbative term, the large corrections in 
the numerator and denominator cancel out with each other (see Eq.(10)). In fact, the 
£(s,m c ) are complicated functions of the energy s and the mass m c , such a crude esti- 
mation maybe underestimate the 0(a s ) corrections, the uncertainties originate from the 
0(a s ) corrections maybe larger. 

In this article, we also neglect the contributions from the perturbative corrections 
0(a™). Those perturbative corrections can be taken into account in the leading logarithmic 
approximations through anomalous dimension factors. After the Borel transform, the 
effects of those corrections are to multiply each term on the operator product expansion 
side by the factor, 

r2M 2r J- r o n 

(12) 



a s (M 2 



a s {p?) 

where the Tj is the anomalous dimension of the scalar interpolating current J(x), the To n 
is the anomalous dimension of the local operator O n (0) in the operator product expansion, 

t{j(x)J\0)} = C n (x)O n (0), (13) 

here the C n (x) is the corresponding Wilson coefficient. 

We carry out the operator product expansion at a special energy scale, say /x = 1 GeV, 
and can not smear the scale dependence by evolving the operator product expansion side 
to the energy scale M through Eq.(12) as the anomalous dimension of the scalar current 
J(x) is unknown. Furthermore, the anomalous dimensions of the high dimensional local 
operators have not been calculated yet, and their values are poorly known. In this article, 



we set the factor 



q4m2) 1 2Fj 

as(A» 2 ) 



~ 1, such an approximation maybe result in some scale 
dependence and weaken the prediction ability; further studies are stilled needed. 

In the QCD sum rules, the high dimension vacuum condensates are always factorized 
to lower condensates with vacuum saturation, factorization works well in large N c limit. 
In the real world, N c = 3, there are deviations from the factorable formula, we introduce 
a factor k to parameterize the deviations, 

(ss) 2 , (ss)(sg s aGs) , (sg s aGs) 2 — > k(ss) 2 , K(ss)(sg s aGs) , n(sg s aGs) 2 . (14) 
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In Fig. 6, we show the mass My with variation of the parameter k at the interval k = 
— 2. From the figure, we can see that the value of the My changes quickly at the region 
M 2 < 2.6 GeV 2 , and increase with the k monotonously at the region M 2 < 3.2 GeV 2 . At 
the interval M 2 = (2.6 — 3.0) GeV 2 , the value k = 1 ± 1 leads to an uncertainty about 
50MeV, which is too small to smear the discrepancy between the present prediction and 
the experimental data. In the limit k = 0, which corresponds to neglecting the vacuum 
condensates (ss) 2 , (ss)(sg s aGs) and (sg s aGs) 2 , we obtain the smallest value. It is not 
unexpected. From Fig.2E-2F, we can see that there are cancelations among the vacuum 
condensates (ss) 2 , (ss)(sg s aGs) and (sg s aGs) 2 , the net contribution is rather small. If 
we assume the k has the typical uncertainty of the QCD sum rules, say about 30%, 
the correction is rather mild, we can neglect the uncertainty safely and take k = 1, i.e. 
the factorization works well. In the QCD sum rules for the masses of the p meson and 
nucleon, k > 1 [Hj. If the same value holds for the tetraquark states, the deviation from 
the factorable formula means even larger discrepancy between the present prediction and 
the experimental data. 

The c-quark mass appearing in the perturbative terms (see e.g. Eq.(17)) is usually 
taken to be the pole mass in the QCD sum rules, while the choice of the m c in the leading- 
order coefficients of the higher-dimensional terms is arbitrary [15] . The MS mass m c (m 2 ) 
relates with the pole mass rh through the relation 



where K depends on the flavor number n/. In this article, we take the approximation 
m c ~ rh without the a s corrections for consistency. The value listed in the Particle 
Data Group is rn c (m 2 ) = 1.27^0 ^ GeV [16], it is reasonable to take the value m c = 
m c (l GeV 2 ) = (1.35 ± 0.10) GeV. In Fig. 4, we also present the result with smaller value 
m c = 1.3 GeV, which can move down the central value about 0.06 GeV. The central value 
My = 4.37 GeV is still larger than the D*D* threshold about 150 MeV. 

The QCD sum rules is just a QCD-inspired model, we calculate the ground state mass 
by imposing the two criteria (pole dominance and convergence of the operator product 
expansion) of the QCD sum rules. In fact, we can take smaller threshold parameter sq 
and larger Borel parameter M 2 to reproduce the experimental value by releasing the pole 
dominance condition. 

We usually consult the experimental data in choosing the Borel parameter M 2 and the 
threshold parameter sq. The present experimental knowledge about the phenomenological 
hadronic spectral densities of the multiquark states (irrespective of the molecule type and 
the diquark-antidiquark type) is rather vague, even existence of the multiquark states is 
not confirmed with confidence. The nonet scalar mesons below 1 GeV (the /o(980) and 
ao (980) especially) are good candidates for the tetraquark states. However, they can't 
satisfy the two criteria of the QCD sum rules, and result in a reasonable Borel window. If 
the perturbative terms have the main contribution (in the conventional QCD sum rules, the 
perturbative terms always have the main contribution), we can approximate the spectral 
density with the perturbative term [11] . then take the pole dominance condition, and 
obtain the approximate relation, 




(15) 
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Figure 3: The mass with variation of the Borel parameter M 2 . The notations a, (3, 7, 
A, r, £ and p correspond to the threshold parameters so = 19 GeV 2 , 20 GeV 2 , 21 GeV 2 , 
22 GeV 2 , 23 GeV 2 , 24 GeV 2 and 25 GeV 2 , respectively. 

If the Borel parameter has the typical value M 2 = lGeV 2 , then sq > 4.7 GeV 2 , the 
threshold parameter is too large for the light tetraquark state candidates /o(980), ao(980), 
etc. 

Once the main Fock states of the nonet scalar mesons below 1 GeV 2 are proved to be 
tetraquark states, we can draw the conclusion that the QCD sum rules are not applicable 
for the light tetraquark states. We can either reject the QCD sum rules for the multiquark 
states or release one of the two criteria (pole dominance and convergence of the operator 
product expansion) [8], for example, we can cut the threshold parameters for the /o(980) 
and ao(980) slightly larger than 1 GeV 2 by hand. 

4 Conclusion 

In this article, we assume that there exists a scalar D*D* molecular state in the J/ip(f> 
invariant mass distribution, and study its mass using the QCD sum rules. Our predic- 
tions depend heavily on the two criteria (pole dominance and convergence of the operator 
product expansion) of the QCD sum rules. The numerical result indicates that the mass 
is about My = (4.43 db 0.16) GeV, which is inconsistent with the experimental data. The 
D*D* is probably a virtual state and not related to the meson y(4140). Other possibility, 
such as a hybrid charmonium is not excluded; more experimental data are still needed to 
identify it. 
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Figure 4: The mass with variation of the Borel parameter M 2 . 
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Figure 5: The pole residue with variation of the Borel parameter M 2 . 
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Figure 6: The mass with variation of the parameters k and M 2 , other parameters are 
taken to be the central values. 

Appendix 

The spectral densities at the level of the quark-gluon degrees of freedom: 

Po(s) = — e / da d(3a(3(l — a — (3) 3 (s — m 2 ) 2 (7s 2 — 6sm 2 + fn^) 

10247r b J a . J p . 

+ 10^6 f 1 da I a d(3a(3(l-a-P) 2 (s-m 2 c f(3s-m 2 c ) 
Srr? m f' a f rl— a 

J da J df3{a + m ~ a " p)2[s ~ ™' )2(5s " 2 ™ 2) ' (17) 
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Q m f a f pi.— a 

P(ss)(s) = "^V / da dPaf3(l-a-(3)(Ws 2 -l2sfh 2 c +3mi) 

+^0 L £ ^ d(3af3(s - m 2 c ){2s - m 2 .) 

_ ^{sgsaGs) £ da J 1 - [6(2g _ ^ + s2g{s _ ^ 2)] 



3m c (ss) r a f r 1 -" 
32vr 4 



da I d(i{a + /3)(1 - a - [3){s - m 2 c )(2s - m\) 

at J 0i 



3m c (sg s aGs) f a f f 1 ^ 
128~^ / / dp(a + (3)(3s -2m 



I) 

ft 



m s (sg s aGs) f af . , w ~2 
64^ J dQQ (l-«)(3s-H) 

3m s m 2 c (sg s aGs) [ a f 
+ 32^ I **> (18) 



m r . 



(ss) 2 f a f m 2 c {sg s aGs) 2 f a f ~4 



4vr 2 



da + 



m 2 (ss){sg s aGs) f a f 



8vr 2 

m s m c (ss} 2 f a f 
16vr 2 



da 



64tt 2 M 6 .,„ 

1 H - 

M 2 



da 



dam c 5(s — m c 



5(s — m c 



2 + s5(s — m c 



5m s m c (ss)(sg s aGs) f af 



96vr 2 



f 



da 



i , m c I TO C 

M 2 2M 4 



<S(*-m c ), (19) 
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P{GG) ( s ) 



256vr 4 



r f r 1 
/ da 

JOLi J fii 



' a ^f^+4)( i -«-/ ? ) 3 



(5 1 a 2 



2s — m., H 5(s — m„) 

6 



3m s m c in';. /'" ' , I ' " , , f ^ \ . , j 2 



512vr 4 



m s m\ 

' 5127T 4 



3 /-a 



1-a 



'/< « / d/3 + ) (a + /3)(1 - a - /?) 2 [2 + s5(s - m 2 )] 

l-a 



] rl—a 

j/ da d(3{a + 0){l - a - (3) 2 {lQs 2 - 12sm 2 c + 3fh^) 

512vr 4 ,/ a . y a . 

d/3(a + - a - /3)(s - m 2 )(2s - m 2 ) 



'Oli •> Bi 

^ pctf pi— a 



256vr 4 



da 

J Pi 

a 



1m m r a f rl—a 

SmsTnc 1 da I d(3(l - a - P)(3s - 2m 2 ) 



1287T 4 j ai JPi 
m s ml { Ss) n da ( l - a di 

hi 

5{s-m 2 c ) 



96vr 2 



~ 2 ~ 4 

1 + tt^ + 



M 2 2M 4 
m s m 2 (ss) f a f 

1927T 2 



my rl—a 

/ da / 



1 + 



m c 

M 2 



5(s — m 



m s m 2 c (sg s aGs) f a f f 1 a ( a \ ~ 2 n 



m s mf(ss) f af 



48tt 2 M 2 
m^.(ss) 



/ da / 



a 4 Jft 

1 + ^ 1^-^) 



l-a 



1927T 2 

~ 2 
m r 

1 H - 

M 2 



/ da/ d/?( — + ^ ) (a + /3)(l-a-/3) 



a 3 /3 3 



5(s — m 



(3)m 2 5(s — m 



ml(sg s aGs) f a f 
768vr 2 M 4 

- r * /r * 3 (? + ^' <i_a " ,j) ' 2+ si<s - ™ 3] 



m c (sg s aGs) 



d P ( — + 4 

A 

l-a 



/3 2 a 2 



/•a/ 

/ da 

Jai Jp 



256vr 2 
m s m 2 (ss) 

167T 2 

+!?^!£l da T " d/3 [ 2 + - m 2 )] , 



1 + 



M 2 

2\ 



— m 



32vr 2 



Pi 



(20) 



12 



P(GG}( S ) 



m 4 (ss)' 
72M 4 



da 



1 

s + 



a 3 (1 — a) 3 



m s mf(sg s o-Gs} f af 



+ 



192vr 2 M 4 
m s m 2 c {sg s (jGs ) [ a f 
1152vr 2 M 4 



da 



+ 



da 



+ 



m s m 3 (ss)' 
288M 4 

m 2 (ss) 2 
2AM 2 



m s m 2 c (sg s oG s) [ af 

64?r 2 M 2 
m s m c (ss) 2 



da 

da 

' 1 

a'- 



1 

-7? + 



<5(s -m c ) 
1 



a 3 (1 — a) c * 



1 — a 



+ 



a 



a" 



+ 



1 



1 

— + 

a 



(1-a) 3 
1 



1-a) 2 
1 - 



«5(s - m c ) 



m c 5(s-m c ) 



M 2 



5(s-m c ) 



2 ■ (i 



da 



1 



a 



+ 



5{s-m c ) 
1 



96M 4 
m s {ss) 



da 



a 

1 - a 



+ 



(1-a) 2 
a 



384vr 2 

m c (sg s aGs 
128vr 2 



da 



a' (1 — a) 2 

-2 



J(s - m c ) 



~2 ~2 

m c 5(s — m c ) 



2 + s5(s — m c ) 

~2 



where a/ 



i+Vi 



1-V1- 



da 



A 



i + 



M 2 



<5(a -m c ). 



: ~ 9 (a+/3)mi ~ 



a/3 ' "*c — a(l-a) ' 



(21) 
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